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ABSTRACT 

The migration of the giant planets due to the scattering of planetesimals causes powerful resonances to move through the asteroid 
belt and the terrestrial planet region. Exactly when and how the giant planets migrated is not well known. In this paper we present 
results of an investigation of the formation of the terrestrial planets during and after the migration of the giant planets. The latter is 
assumed to have occurred immediately after the dissipation of the nebular disk - i.e. "early" with respect to the timing of the Late 
Heavy Bombardment (LHB). The presumed cause of our modeled early migration of the giant planets is angular mometum transfer 
between the planets and scattered planetesimals. 

Our model forms the terrestrial planets from a disk of material which stretchs from 0.3-4.0 AU, evenly split in mass between planetes- 
imals and planetary embryos. Jupiter and Saturn are initially at 5.4 and 8.7 AU respectively, on orbits with eccentricities comparable 
to the current ones, and migrate to 5.2 and 9.4 AU with an e-folding time of 5 Myr. 

Unfortunately, the terrestrial planets formed in the simulations are not good analogs for the current solar system, with Mars typically 
being much too massive. Moreover, the final distribution of the planetesimals remaining in the asteroid belt is inconsistent with the 
observed distribution of asteroids. This argues that, even if giant planet migration had occurred early, the real evolution of the giant 
planets would have to have been of the "jumping- Jupiter" type, i.e. the increase in orbital separation between Jupiter and Saturn 
had to be dominated by encounters between Jupiter and a third, Neptune-mass planet. This result was already demonstrated for late 
migrations occuring at the LHB time by Brasser et al. (2009) and Morbidelli et al. (2010), and this paper shows their conclusions hold 
for early migration as well. 

Key words, planet formation - asteroid belt 



1. Introduction 

The formation of the terrestrial planets is expected to have oc- 
curred from a disk of planetesimals in two steps. In the first 
step, Moon to Mars-size "planetary embryos" formed by run- 
away and oligarchic accretion (Greenberg et al. 1978; Wetherill 
& Stewart 1993; Kokubo & Ida 1998). In the second step, the 
terrestrial planets formed by high-velocity collisions among the 
planetary embryos (Chambers & Wetherill 1998; Agnor et al. 
1999; Chambers 2001; Raymond et al. 2004, 2005, 2006, 2007; 
O'Brien et al. 2006; Kenyon & Bromley 2006). 

The most comprehensive effort to date in modeling terrestrial 
planet formation (Raymond et al. 2009) focused on 5 constraints 
of the terrestrial planets: 1. the orbits, particularly the small ec- 
centricities, 2. the masses, with the small mass of Mars the most 
difficult to match, 3. formation timescales, 4. bulk structure of 
the asteroid belt and 5. the water content of Earth. Despite suc- 
cess with some of these constraints in each simulation, no sim- 
ulation satisfied all the constraints simultaneously. For the sim- 
ulations with fully formed Jupiter and Saturn on nearly circu- 
lar orbits, the constraint consistently missed is the small mass 
of Mars. A Mars of the correct size is only obtained in simu- 
lations where the giant planets are on orbits with current semi- 
major axes but much larger eccentricities. This scenario, how- 
ever, raises the problem of not allowing any water delivery to 
Earth from material in the outer asteroid belt region. The size 
of Mars has been a consistent problem for previous works with 
giant planets assumed on current orbits and disks of planetesi- 
mals and embryos stretching from ~0. 5^4.0 AU (Chambers & 



Wetherill 1998; O'Brien et al. 2006), or even only up to 1.5 or 
2.0 AU (Kokubo et al. 2006; Chambers 2001). 

However, Hansen (2009) had great success creating analogs 
of Mars in simulations which begin with a narrow annulus of 
planetary embryos between 0.7 and 1.0 AU. In these simulations 
both Mercury and Mars are formed from material that is scat- 
tered out of the original annulus by the growing Earth and Venus 
analogs. In addition, the orbits of the Earth and Venus analogs 
have eccentricities and inclinations similar to those observed to- 
day and the accretion timescales are in agreement, although on 
the low side, with the ages of the Earth-Moon system deduced 
from the 182 Hf - 182 W chronometer. This model points to the 
need for a truncated planetesimal disk at, or near, the beginning 
of the process of terrestrial planet formation. The origin of this 
truncation remains to be understood. Similarly, it remains to be 
clarified how the truncation of the disk of planetesimals at 1 AU 
can be compatible with the existence of asteroids in the 2-4 AU 
region. 

Nagasawa et al. (2005) and Thommes et al. (2008) effec- 
tively produced a cut in the planetesimal distribution at 1 .5 AU 
by assuming that the giant planets were originally on their cur- 
rent orbits and that secular resonances swept through the aster- 
oid belt during gas-dissipation. However, the assumption that 
the giant planets orbits had their current semimajor axes when 
the gas was still present is no longer supported. When embed- 
ded in a gas disk, planets migrate relative to each other until a 
resonance configuration is achieved (Peale & Lee 2002; Kley et 
al. 2009; Ferraz-Mello et al. 2003; Masset & Snellgrove 2001; 
Morbidelli & Crida 2007; Pierens & Nelson 2008). Thus it is 
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believed that the giant planets were in resonance with each other 
when the gas disk disappeared (Morbidelli et al. 2007; Thommes 
et al. 2008; Batygin & Brown 2010) which causes problems in 
understanding the consequences of the Thommes et al. (2008) 
model. Moreover, the Nagasawa et al. (2005) and Thommes et 
al. (2008) simulations produce the terrestrial planets too quickly 
(~ 10 Myr), compared to the timing of moon formation indi- 
cated by the 182 Hf - 182 W chronometer (> 30 Myr and most 
likely > 50 Myr; Kleine et al. 2009) and they completely deplete 
the asteroid belt by the combination of resonance sweeping and 
gas-drag (see also Morishima et al. 2010, for a discussion). 

The resonant configuration of the planets in a gas disk is 
extremely different from the orbital configuration observed to- 
day. Planetesimal-driven migration is believed to be the mech- 
anism by which the giant planets acquired their current orbits 
after the gas-disk dissipation. In fact, work by Fernandez & 
Ip (1984) found that Uranus and Neptune have to migrate out- 
ward through the exchange of angular momentum with planetes- 
imals that, largely, they scatter inward. Similarly, Saturn suffers 
the same fate of outward migration, though Jupiter migrates in- 
ward as it ejects the planetesimals from the solar system. The 
timescale for planetesimal-driven migration of the giant planets 
depends on the distribution of the planetesimals in the planet- 
crossing region. It is typically 10 My, with 5 My as the lower 
bound (Morbidelli et al. 2010). Close encounters between pairs 
of giant planets might also have contributed in increasing the or- 
bital separations among the giant planets themselves (Thommes 
et al. 1999; Tsiganis et al. 2005; Morbidelli et al. 2007; Brasser 
et al. 2009; Batygin & Brown 2010). Beyond the consequences 
for the scattered planetesimals, the migration of the giant plan- 
ets affects the evolution of the solar system on a much larger 
scale, through the sweeping of planetary resonances through the 
asteroid belt region. 

The chronology of giant planet migration is important for the 
general evolution of the solar system, including the formation of 
the terrestrial planets. It has been recently proposed (Levison et 
al. 2001; Gomes et al. 2005; Strom et al. 2005) that the migra- 
tion of the giant planets is directly linked in time with the so- 
called "Late Heavy Bombardment" (LHB) of the terrestrial plan- 
ets (Tera et al. 1974; Ryder 2000, 2002; Kring & Cohen 2002). 
If this is true, then the migration of the giant planets should have 
occurred well after the formation of the terrestrial planets. In 
fact, the radioactive chronometers show that the terrestrial plan- 
ets were completely formed 100 Myr after the condensation of 
the oldest solids of the solar system (the so-called calcium al- 
luminum inclusions, which solidified 4.568 Gyr ago; Bouvier 
et al. 2007; Burkhardt et al. 2008), whereas the LHB occurred 
3.9-3.8 Gyr ago. Thus the terrestrial planets should have formed 
when the giant planets were still on their pre-LHB orbits: reso- 
nant and quasi-circular. However, the simulations of Raymond 
et al. (2009) fail to produce good terrestrial planet analogs when 
using these pre-LHB orbits. 

The alternative possibility is that giant planet-migration oc- 
curred as soon as the gas-disk disappeared. In this case, it can- 
not be a cause of the LHB (and an alternative explanation for 
the LHB needs to be found; see for instance Chambers 2007). 
However, in this case giant planet migration would occur while 
the terrestrial planets are forming, and this could change the out- 
come of the terrestrial planet formation process. In particular, 
it is well known that, as Jupiter and Saturn migrate, the strong 
vg secular resonance sweeps through the asteroid belt down to 
~ 2 AU (Gomes 1997). The V6 resonance occurs when the pre- 
cession rate of the longitude of perihelion of the orbit of an as- 
teroid is equal to the mean precession rate of the longitude of 



perihelion of Saturn, and it affects the asteroids' eccentrcities. If 
the giant planet migration occurs on a timescale of 5-10 Myr, 
typical of planetesimal-driven migration, then the ve resonance 
severely depletes the asteroid belt region (Levison et al. 2001; 
Morbidelli et al. 2010). This can effectively truncate the disk of 
planetesimals and planetary embryos, leaving it with an outer 
edge at about 1 .5 AU. Although the location of this edge is not 
as close to the sun as assumed in Hansen (2009) (1 AU), it might 
nevertheless help in forming a Mars analog, i.e. signficantly less 
massive than the Earth. 

An equally important constraint is the resulting orbital distri- 
bution of planetesimals in the asteroid belt region, between 2-4 
AU. After that region has been depleted of planetesimals and 
embryos by the sweeping resonances, what remains will survive 
without major alteration and should compare favorably with to- 
days large asteroids. Studies of late giant planet migration start 
with an excited asteroid belt, where inclinations already vary 
from 0-20° (Morbidelli et al. 2010), and cannot match the in- 
clination distribution of the inner asteroid belt with 5 Myr or 
longer migration timescales. The early migration presented here 
is different because it occurs immediately after the dissipation 
of the gas disk so that the planetesimal orbits are dynamically 
cold, with inclincations less than 1°. Thus, in principle, an early 
giant planet migration could lead to a different result. Also, the 
embryos will be present, another difference with late migration 
scenarios. 

The purpose of this paper is to investigate, for the first time, 
the effect that an early migration of the giant planets could have 
had on the formation of the terrestrial planets and on the final 
structure of the asteroid belt. In Section 2 we discuss our meth- 
ods and in Section 3 we present our results. The conclusions and 
a discussion on the current state of our understanding of terres- 
trial planet formation will follow in Section 4. 



2. Methods 

We assume in our simulations that the nebular gas has dissipated, 
Jupiter and Saturn have fully formed; in the terrestrial planet and 
asteroid belt region, in the range 0.5-4.0 AU, the planetesimal 
disk has already formed planetary embyros accounting for half 
of its total mass. The lifetime of the circumstellar gas disk is 
observed to be 3-6 Myr, and both Jupiter and Saturn are ex- 
pected to be fully formed by this time (Haisch et al. 2001). The 
timescales for oligarchic growth is similar, with lunar to Mars 
sized embryos growing on million year timescales (Kokubo & 
Ida 1998,2000). 

The numerical simulations are done using SyMBA, a sym- 
plectic N-body integrator modified to handle close encounters 
(Duncan et al. 1998). In our model, the planetary embryos in- 
teract with each other; the planetesimals interact with the em- 
bryos but not with themselves; all particles interact with the gi- 
ant planets and, except when specified (explained further below), 
the giant planets feel the gravity of embryos and planetesimals. 
Collisions between two bodies result in a merger conserving lin- 
ear momentum. It has been demonstrated by Kokubo & Genda 
(2010) that this a priori assumption of simple accretion does not 
significantly affect the results. The SyMBA code has already 
been used extensively in terrestrial planet formation simulations 
(Agnor et al. 1999; Levison & Agnor 2003; O'Brien et al. 2006; 
McNeil et al. 2005). 
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2.1. Protoplanetary disk 

The initial protoplanetary disks are taken directly from O'Brien 
et al. (2006), which themselves were based on those of 
Chambers (2001). The O'Brien et al. (2006) study produced 
some of the best matches for terrestrial planets and by using 
similar intial conditions allows a direct comparison. The initial 
conditions are based on a "minimum mass" solar nebula, with a 
steep surface denstiy profile. The solid mass is shared between 
many small planetesimals and a small number of large bodies, 
the embryos, as suggested by runaway/oligarchic growth simu- 
lations (Kokubo & Ida 1998, Kokubo & Ida 2000). In theory, 
it is possible that, by the time the gas disappears from the disk 
(which cooresponds to time zero in our simulations) the plane- 
tary embryos in the terrestrial planet region could have grown 
larger than the mass of Mars. However, the current mass of Mars 
seemingly excludes this possibility, and argues for masses to 
have been martian or sub-martian in mass. 

The surface density profile is Z(r) = So(j^j) -3 ^ 2 , where 
So = 8 g crrT 2 . The distribution of material drops linearly be- 
tween 0.7 and 0.3 AU. Half of the mass is in the large bodies, 
of which there were either 25 embryos, each of 0.0933 Earth 
masses (M @ ) or 50 embryos of 0.0467 M m . The small bodies 
are 1/40 as massive as the large embyros, or 1/20 as massive as 
the small embryos. For all test cases the embryos are spaced be- 
tween 4-10 mutual hill radii at the beginning of the simulations. 
In some tests, the smaller planetesimals with an initial semimajor 
axis larger than 2.0 were cloned into two particles with identi- 
cal semi-major axis, half the mass in each, and different random 
eccentricities and inclinations (noted as 'Double Asteroids' in 
Table 1.). The initial eccentricities and inclinations were selected 
randomly in the range of 0-0.01 and 0-0.5 degrees respectively. 
Thus the initial mass of the disk consisted of 2.6 M e located 
inside of 2 AU and a total mass of 4.7 M m . 

2.2. Giant planets and migration 

In all tests Jupiter and Saturn were started on orbits closer to 
each other than at the present time, i.e. with semimajor axes of 
5.4 and 8.7 AU, respectively. These initial orbits are just be- 
yond their mutual 1:2 mean motion resonance, i.e. the corrse- 
ponding ratio of orbital periods of Saturn and Jupiter is slightly 
larger than 2. Even if the giant planets should have started from 
a resonant configuration - probably the 2:3 resonance (Masset 
& Snellgrove 2001; Pierens & Nelson 2008) - it is known that 
secular resonance sweeping through the asteroid belt is impor- 
tant only when the planets' orbital period ratio is larger than 2 
(Gomes 1997; Brasser et al. 2009). Thus, our choice of the ini- 
tial orbits of Jupiter and Saturn is appropriate for the purposes 
of this study. 

Each planet was forced to migrate by imposing a change to 
their orbital velocities that evolves with time t as: 

v{t) = vo + Av[l - exp(-f/r)] 

appropriate Av to achieve the required change in semimajor axis, 
and t — 5 Myr. The latter is the minimum timescale at which 
planetesimal-driven migration can occur, simply due to the life- 
time of planetesimals in the giant planet crossing region, as dis- 
cussed extensively in Morbidelli et al. (2010). Longer timescales 
are possible, but previous work has shown that fast timescales 
affect the asteroid belt region less, and since terrestrial planet 
formation timescales are in the 10's of millions of years, more 
rapid migration has a greater chance of affecting the accretion 
of Mars. Thus, we think that restricting ourselves to the 5 My 
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Fig. 1. Example of idealized migration for a system with only 
Jupiter and Saturn, ending with orbits very close to the current 
ones. Panel (a) shows the semimajor axis of Jupiter, (b) the ec- 
centricity of Jupiter, (c) the semimajor axis of Saturn and (d) the 
eccentricity of Saturn, all plotted as a function of time in years. 



timescale is sufficient, as this timescale is the most favorable for 
these purposes. 

If the motion of the giant planets was not affected by the 
other bodies in the system, the evolution of the eccentricities and 
inclinations would not change much during migration (Brasser 
et al., 2009, and Fig. 1). Thus it is relatively simple to find initial 
conditions that lead to the final orbital configurations with ec- 
centricities and inclinations with mean values and amplitude of 
oscillations similar to current one. In fact, as shown in Brasser 
et al. (2009), the initial values (e y ,e s ) = (0.012,0.035) and 
(ij,is) = (0.23°, 1.19°), after migration, lead to eccentricities 
and inclinations whose mean values and amplitudes of oscil- 
lation closely resemble those characterizing the current secular 
dynamics of the giant planets (see Fig.[T|). 

In our case, however, as the giant planets migrate, they scat- 
ter planetesimals and planetary embryos, and their orbits are af- 
fected in response. Thus, the final orbits are not exactly like those 
of Fig. [1] Typically, for instance, the eccentricities and inclina- 
tions of the planets are damped, and their relative migration is 
slightly more pronounced than it was intended to be. Thus, we 
tried to modify the initial eccentricities of Jupiter and Saturn and 
the values of Av in order to achieve final orbits as similar as pos- 
sible to those of Fig. Q] However, while the effect of planetesi- 
mals on the planets is statistically the same from simulation to 
simulation, (and so can be accounted for by modifying the initial 
conditions of the planets), the effects of embryos are dominated 
by single stochastic events. Thus, it is not possible to find plan- 
etary initial conditions that lead systematically to good final or- 
bits. In some cases the final orbits are reasonably close to those 
of the current system, but in many cases they are not. In total we 
performed 30 simulations. We discarded the simulations with 
unsuccessful final orbits, and kept only those (9/30) that lead 
to orbits resembling the current ones. These successful runs are 
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called hereafter "normal migration simulations". Our criterion 
for discriminating good from bad final orbits was determined af- 
ter the 15 Myr of migration, and the semimajor axis, eccentricity 
and oscillation in eccentricity (Ae) were the factors examined. 

Jupiter's orbit must have had \a—aj\ < 0.05, \e—ej\ < 0.0156, and 1 " ^g>^ 

\Ae-Aej\ < 0.0164, while Saturn's orbit required \a-a s \ < 0.075, ^ Ifca^* n 

\e-e s \ < 0.0252, and Ae-Ae s | < 0.0256. ^ ~ ^ 

ag & 

We complemented our normal migration simulations with ^ ^ a a 

what we call hereafter 'perfect migration' cases. In these simu- 1 
lations, the planetesimals and embryos do not have any direct ef- a "it— 
feet on the giant planets, even during close encounters. However, 1 

their indirect effects cannot be suppressed (specifically the H sun " ^ — b 

term from eq. 32b. in Duncan et al. 1998) , but in principle they o.i - ■ A = 

are weaker. Thus the migration of the giant planets, starting with 

the initial conditions from Brasser et al. 2009 (as in Fig. [TJ, met — q- 
the above criteria in 3 out of 4 simulations. The giant planets 
had the full gravitational affect on the planetesimals and em- 
bryos throughout these simulations, and the mutual effects be- 
tween planetesimals and embryos remained unchanged. 



1.5 

semimajor axis {AU) 



3. Results 

We present the results of 12 simulations of terrestrial planet for- 
mation each covering 150 Myr. Of these runs, 9 are normal mi- 
gration simulations and 3 'perfect migration' simulations (all 
simulations are listed in Table 1. and refered to by run name, 
"Test31" etc., throughout). These two sets of simulations had 
qualitative and quantitative similarities and are thus discussed 
at the same time and combined in the figures. First, the result- 
ing planets are compared with the current terrestrial planets, fol- 
lowed by a look at the consequences the migration has on the 
structure of the asteroid belt. 



3.1. The planets 

Results for these simulations are summarized in Fig. |2l where 
the final masses and semimajor axes of our synthetic planets are 
compared to those of the real terrestrial planets (see also Table 
[TJ. The trend is similar to that found in previous works (see for 
instance Chambers et al. 2001), where the masses and locations 
of Earth and Venus are nearly matched by a number of different 
simulations, but most planets just exterior to Earth, near a ~1.5 
AU are at least 3 times more massive than Mars. However, a 
handful of planets close to 1.8 AU were of similar mass to Mars. 
Of note, Test31 had two ~Mars-mass bodies, at 1.2 and 2.4 AU, 
with an Earth mass planet at 1.52 AU. Test54, the only one of 
four simulations starting with the smaller embryos with success- 
ful migration, produced a sub-Mars mass body at 1.89 AU, just 
at the edge of the current day asteroid belt. The Ran4 simulation 
produced a body within 50% of Mars' mass at 1.71 AU, though 
it had a high eccentricity above 0.13 and was a member of a 3 
planet system. In general, planets produced at around 1.5 AU 
were ~ 5 times more massive than Mars, and Mars-mass bodies 
were typically only found beyond 1.7 AU. 

The total number of planets produced in each simulation is 
not systematically consistent with the real terrestrial planet sys- 
tem. Only two simulations produce 4 planets, where we define 
a "planet" as any embryo-sized or larger body with a semimajor 
axis less than 2.0 AU. Most simulations had 3 planets at the end, 
while one produced 5 planets. A common metric for measur- 



Fig. 2. The final mass (M m ) for each planet produced in our sim- 
ulations is plotted as a function of the planet's semimajor axis. 
The horizontal error bars show the locations of the perihelion 
and aphelion of the cooresponding orbit. The open squares re- 
fer to the planets produced in the normal migration simulations, 
the open circles to the planets produced in the run with twice as 
many half-sized embryos, and the open triangles to those pro- 
duced in the 'perfect migration' simulations; the solid squares 
represent the real terrestrial planets. 

ing the distribution of mass among multiple planets is the radial 
mass concentration statistic (RMC), defined as 



RMC = max] 



ZMj\log 10 (a/aj)y 



(D 



where Mi and aj are the mass and semimajor axis of planet j 
(Chambers, 2001). The bracketed function is calculated for dif- 
ferent a in the region where the terrestrial planets form. The 
RMC is infinite for a single planet system, and decreases as mass 
is spread among multiple planets over a range of semimajor axes. 
The current value of RMC for the solar system is 89.9. For all 
but one simulation the RMC value is below the current solar sys- 
tems value, largely due to the large mass concentrated in a Mars- 
analog orbit (we did not include the two embryos stranded in the 
asteroid belt region in these calculations, one in Test3 1 and one 
in TestPM24). The single simulation with a larger RMC value 
did not have a Mars analog, and thus the mass was contained in 
a smaller semimajor axis range. 

The terrestrial planets have low eccentricities and inclina- 
tions, Earth and Venus both have e < 0.02 and i < 3°, proper- 
ties which has proved difficult to match in accretion simulations. 
O'Brien et al. (2006) and Morishima et al. (2008) reproduced 
low eccentricities and inclinations largely due to remaining plan- 
etesimals which damp the orbital excitation of the planets. A 
metric used as a diagnostic of the degree of success of the simu- 
lations in reproducing the dynamical excitation of the terrestrial 
planets is the angular momentum deficit (AMD; Laskar 1997): 



AMD 



ZjMj^j(l-cos(ij)JT^efj 



(2) 
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Fig. 3. Evolution of the system over time, showing the clearing 
of the asteroid belt region with inclination plotted as a function 
of semimajor axis. The open boxes are planetesimals on orbits 
within the current asteroid belt region, the crosses are planetes- 
imals elsewhere, and the open circles are embryos or planets 
scaled in relation to their diameters. The simulation is Test31. 



where Mi and aj are again the mass and semimajor axis and 
ij and ej are the inclination and eccentricity of planet j. The 
AMD of the current solar system is 0.0014. The AMD for our 
simulations ranged from 0.001 1 to 0.01 13. The plantesimal disk 
used in these simulation is based on that from O'Brien (2006), 
and is therefore not surprising that some AMDs are consistent 
with the solar system value. Simulation PM22 is the one with 
the largest final AMD, because it produced an Earth-analog with 
a 10° inclination. 

Figures[3]and|4]show snapshots of two systems evolving over 
time. Of interest is the radial clearing caused by the movement 
of the giant planets and the sweeping of their resonances, partic- 
ularly the Vft resonance. This clearing progresses from the outer 
edge of the disk towards the sun, following the migration of the 
V6 resonance, and stops at ~ 2 AU, which is the final location of 
this resonance when the giant planets reach their current orbits. 
Thus, the region of a > 2.0 is almost entirely cleared of material 
in 10 Myr, with only handfuls of planetesimals surviving and a 
single embryo. At 3 Myr, only a > 2.5 AU is largely cleared. 

As seen in Figure [5] the accretion of embryos for the Mars 
analogs (where a Mars analog is defined as the largest body be- 
tween 1 .2-2.0 AU) begins immediately with ~2 Mars-mass typ- 
ically being reached in only 2 Myr (note that Figure [5] shows 12 
growth curves, as there are two planets displayed for Test54). 
Nine of the 1 1 Mars analogs have reached 0.2 M e by 3 Myr. 
At 10 Myr 6 of the 1 1 have reached 0.3 M®, and by 30 Myr 10 
of 1 1 are above 0.3 M m , or ~3 MMars- One might wonder if our 
inability to produce a small Mars analog is due to the fact that, 
in all but one of the presented simulations (Test54 is the excep- 
tion), the planetary embryos are initially ~ one Mars mass. This 
is not regarded as a problem for the following reasons. First, the 
Mars analogs with semimajor axes near that of Mars, near 1.5 
AU, typically accreted 4 or 5 embryos; thus they consistently 
accreted much more mass than Mars, and are not simply the 



Fig. 4. Same as Fig. [3] but for simulation Test54, which started 
from 50 embryos of 0.0467 M e instead of 25 embryos twice as 
massive. 



result of a chance accretion between two Mars mass embryos. 
Second, only two of the 1 1 Mars analogs did not accrete another 
embryo, in Test54 and Ran4, but both had semimajor axes larger 
than 1.7 AU, well beyond the current orbit of Mars. Third, our 
single successful normal migration simulation that started with 
half-Mars mass embryos also produced an Earth mass planet at 
1.2 AU. This planet was already two-mars masses in 5 Myr (no- 
tice that in the same simulation one embryo escaped all colli- 
sions with other embryos and therefore remained well below the 
mass of Mars - see Figure|5]- but this object ended up at 1.9 AU, 
well beyond the real position of Mars). Finally, previous works 
(Chambers 2001; Raymond et al. 2009; Morishima et al. 2010 to 
quote a few) which started with embryos significantly less mas- 
sive than Mars met the same Mars-mass problem found here. 
The similarities between our work and previous in terms of the 
mass distribution of the synthetic planets as a function of semi- 
major axis suggest that the giant planet migration does not af- 
fect significantly the terrestrial planet accretion process. Thus it 
is unlikely that small changes in the adopted evolution pattern 
of the giant planets could lead to significantly different results. 
Therefore the initial conditions do not appear to be at fault for 
the failure to match the mass of Mars. 

The reason for which the Mars analog consistently grows too 
massive is twofold. First, they grow fast (in a few million years, 
as shown in Figure |5}, compared with the timescale required to 
effectively truncate the disk at ~ 2 AU (10 Myr, as shown in 
Figures [3] and @J. Second, the truncation of the disk caused by 
the sweeping of the V(, resonance is not sunward enough: the 
final edge is approximately at 2 AU, whereas an edge at ~ 1 AU 
is needed (Hansen 2009; Kokubo et al. 2006; Chambers 2001). 

Figure [6] shows the final incination vs. semimajor axis distri- 
butions of all our simulations (respectively, 'normal' and 'per- 
fect' ones). The sizes of the symbols representing the planets are 
proportional to the cubic roots of their masses. Again, the prob- 
lem of the mass of Mars stands out. 
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Fig. 5. Mass growth of the Mars analogs for all simulations plot- 
ted as a function of time. The most massive Mars analogs ex- 
ceed the mass of Mars (0. 1 1 M & ) in only 2-3 Myr, and then in 
the next 10-20 Myr continue to grow to their final sizes, ending 
many times more massive than Mars. The two lines starting from 
~ 0.05 M e are for two planets of simulation 'Test54', the only 
successful normal simulation that started with half-Mars mass 
embryos. The bold line shows the mass growth of the planet end- 
ing at ~ 1.2 AU; the thin line the planet ending at ~ 1 .9 AU. 



3.2. The asteroid belt 

In the previous section we have shown that the an early sweeping 
of secular resonances through the asteroid belt is not useful to 
solve the small-Mars problem. Here we address the question of 
other observational constraints. For this purpose, in this section 
we turn to the asteroid belt, whose orbital distribution is very 
sensitive to the effects of resonance sweeping (Gomes 1997; 
Nagasawa et al. 2000; Minton & Malhotra 2009; Morbidelli et 
al. 2010). 

Morbidelli et al. (2010) have shown that the properties of 
the asteroid belt after the slow migration of the giant planets are 
largely incompatible with the current structure of the asteroid 
belt. However, they assumed that the migration of the giant plan- 
ets occurred late, after the completion of the process of terrestrial 
planet accretion and after the primordial depletion/dynamical 
excitation of the asteroid belt. Thus, that work does not exclude 
the possibility of an early migration. In fact, the outcome of an 
early migration could be very different from that of a late migra- 
tion for two reasons: first, the initial orbits of the plantesimals are 
quasi-circular and co-planar in the early migration case whereas 
they are dynamically excited in the late migration case, which is 
an important difference; second, planetary embryos reside in, or 
cross, the asteroid belt region during the early time of terrestrial 
planet formation, and this process has the potential of erasing 
some of the currently unobserved signatures of resonace sweep- 
ing. 

To compare the planetesimal distribution obtained in our 
simulations with the "real" asteroid population, we focus on as- 
teroids larger than ~50 km in diameter, as in previous works 
(Petit et al. 2001; Minton & Malhotra 2009; Morbidelli et al. 




semimajor axis (AL) 



Fig. 6. Endstates of all simulations with the inclination plotted as 
a function of the semimajor axis with asteroids as open squares, 
non-asteroid planetesimals as crosses and embryos/planets as 
open circles scaled by their mass to the 1/3 power. 



2010). These bodies are a reliable tracer of the structure of the 
asteroid belt that resulted from the primordial sculpting pro- 
cesses), as they are too large to have their orbits altered signifi- 
cantly by the thermal Yarkovsky effect or by collisions (see Fig. 
[7]). Moreover, their orbital distribution (see top panel of Fig [7} 
is not affected by observational biases, because all bodies of this 
size are known (Jedicke et al. 2002). 

The final distribution of the planetesimals residing in the as- 
teroid belt in our 12 simulations is shown in the bottom panel 
of Fig [7] As can be seen, the difference in orbital distribution 
between the real belt and that resulting from the giant planet mi- 
gration process is striking. 

A simple metric used in Morbidelli et al. (2010) to quantify 
the difference in orbital distributions between the real and the 
synthetic belts is the ratio of asteroids above and below the lo- 
cation of the ve secular resonance with semimajor axis below 
2.8 AU. The current day value for asteroids with a diameter 
above 50 km is 0.07. Combining together all the surviving plan- 
etesimals from all our 12 simulations results in a 67/13 ratio, 
in stark contrast to the current value. Thus, our result is quali- 
tatively similar to that of Morbidelli et al. (2010), even though 
our resulting ratio is much larger than that obtained in that work 
(close to 1/1). 

The reason for the large ratio obtained in migration simula- 
tions, as discussed in Morbidelli et al. (2010), is that the migra- 
tion of the giant planets forces the v& and vi6 secular resonances 
to move Sun-ward. More precisely, if the orbital separation of 
Jupiter and Saturn increased by more than 1 AU (as predicted 
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Fig. 7. (Top) The inclination of current day asteroids with abso- 
lute magnitude H < 9.7, corresponding to D > 50 km, plotted as 
a function of their semimajor axis. The long-dashed lines show 
the location of the major mean motion resonances with Jupiter 
and the short-dashed curves the location of the V6 and vi6 secular 
resonances. (Bottom) Surviving planetesimals from the 12 sim- 
ulations, showing a strong depletion of low inclination bodies in 
the inner part of the asteroid belt region. 



by all models and enacted in our simulations), the ve resonance 
sweeps the entire asteroid belt as it moves inwards from 4.5 AU 
to 2 AU; meanwhile the vie resonance sweeps the belt inside of 
2.8 AU to its current location at 1.9 AU (Gomes et al. 1997). In 
the inner asteroid belt, the vk, resonance sweeps first and the ve 
resonance sweeps second. The vi6 resonance occurs when the 
precession rate of the longitude of the node of the orbit of an 
asteroid is equal to the precession rate of the node of the orbit 
of Jupiter, and it affects the asteroid's orbital inclination. Given 
the characteristic shape of the V(, resonance location in the (a, i) 
plane (see Fig |7), the asteroids that acquire large enough incli- 
nation when they are swept by the resonance, avoid being 
swept by V(,; thus, their eccentricities are not affected and they 
remain stable. Conversely, the asteroids that remain at low-to- 
moderate inclinations after the vi6 sweeping are then swept by 
the V6 resonance and their eccentricities become large enough 
to start crossing the terrestrial planet region. These bodies are 
ultimately removed by the interaction with the (growing) terres- 
trial planets. This process favors the survival of high-inclination 
asteroids (above the current location of the V6 resonance) over 
low-inclination asteroids and explains the large ratio between 
these two populations obtained in the resonance sweeping sim- 
ulations. This ratio is larger in our simulations than in those of 
Morbidelli et al. (2010), because the initial orbits of planetesi- 
mals and embryos in our case have small inclinations and ec- 
centricities. Consequently, the secular resonance sweeping can 
only increase eccentricities and inclinations. Conversely, in the 
Morbidelli et al. (2010) simulations, the initial orbits covered a 
wide range of eccentricities and inclinations. Large eccentrici- 
ties or inclinations can be decreased by the secular resonance 
sweeping. Thus, more objects could remain at low-to-moderate 



inclinations after the vi6 sweeping and fewer objects were re- 
moved by the v& sweeping than in our case. 

We conclude from our simulations that the migration of the 
giant planets with an e-folding time of 5 Myr (or longer, as the 
effects of secular resonance sweeping increases with increasing 
migration timescale) is inconsistent with the current structure of 
the asteroid belt, even if it occurred early. In fact, our simulations 
provide evidence that the planetary embryos crossing the aster- 
oid belt during the process of formation of the terrestrial plan- 
ets are not able to re-shuffle the asteroid orbital distribution and 
erase the dramatic scars produced by secular resonance sweep- 
ing. 

4. Discussion and Conclusions 

This paper has investigated the effects of early giant planet mi- 
gration on the inner disk of planetesimals and planetary em- 
bryos. In the context of solar system formation, "early" is im- 
mediately following the disappearence of the gas disk, which is 
identified as time-zero in our simulations. The giant planets are 
migrated towards their current orbits with a 5 Myr e-folding time 
which is appropriate if the migration is caused by planetesimals 
scattering. 

We have shown that the sweeping of secular resonances, 
driven by giant planet migration, truncates the mass distribu- 
tion of the inner disk, providing it with an effective outer edge 
at about 2 AU after about 10 Myr. This edge is too far from 
the Sun and forms too late to assist in the formation of a small 
Mars analog. In fact, Chambers (2001) already showed simi- 
lar results starting from a disk of objects with semi-major axes 
0.3 < a < 2.0 AU, the terrestrial planet accretion process leads 
to the formation of planets that are systematically 3-5 times too 
massive at ~ 1 .5 AU. For completeness, we have continued our 
simulations well beyond the migration timescale of the giant 
planets to follow the accretion of planets in the inner solar sys- 
tem, and we have confirmed Chambers (2001) result. 

Hansen (2009) showed that obtaining planets at ~ 1.5 AU 
that have systematically one Mars mass requires that the disk 
of solid material in the inner solar system had an outer edge at 
about 1 AU. The inability of secular resonance sweeping during 
giant planets migration to create such an edge suggests that a 
different mechanism needs to be found. 

Moreover, our study adds to the continuing inability of mod- 
els with a slow migration of the giant planets, t > 5 Myr, to 
leave an asteroid belt with a reasonable inclination distribution. 
Morbidelli et al. (2010) argued that the only possibility for the 
orbits of Jupiter and Saturn to move away from each other on a 
timescale shorter than 1 Myr is that an ice giant planet (presum- 
ably Uranus or Neptune) is first scattered inwards by Saturn and 
is subsequently scattered outwards by Jupiter, so that the two 
giant planets recoil in opposite directions. They dubbed this a 
"jumping-Jupiter" evolution and showed that in this case the fi- 
nal orbital distribution of the asteroid belt is consistent with that 
observed. Again, Morbidelli et al. (2010) worked in the frame- 
work of a "late" displacement of the orbits of the giant planets. 
Our results in this paper suggest that a jumping-Jupiter evolution 
would also be needed in the framework of an "early" displace- 
ment of the orbits of the giant planets. 

At this point, it is interesting to speculate what the effects of 
an "early" jumping-Jupiter evolution would be on the terrestrial 
planet formation process. In essence, an early jumping-Jupiter 
evolution would bring the giant planets to current orbits at a very 
early time. So, the outcome of the terrestrial planet formation 
process would resemble that of the simulations of Raymond et 
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al. (2009) with giant planets initially with their current orbital 
configuration, labelled 'EJS' in that work. In these simulations, 
though, (see their Fig. 10), the Mars analog is, again, systemati- 
cally too big. It is questionable whether a jumping -Jupiter evolu- 
tion could bring the giant planets onto orbits with current semi- 
major axes but larger eccentricities, as required in the most suc- 
cessful simulations of Raymond et al. (2009), labelled 'EEJS'. 
However, even though jumping -Jupiter evolutions satisfying this 
requirement were found, it is important to note that all of the out- 
comes of the EEJS simulations of Raymond et al. (2009). While 
producing a small Mars in several cases, the EEJS simulations 
failed in general to bring enough water to the terrestrial planets, 
formed the Earth too early compared to the nominal timescale of 
50 Myr and left the terrestrial planets on orbits too dynamically 
excited. For all these reasons an early jumping-Jupiter evolution 
is not a promising venue to pursue for a successful model of 
terrestrial planet formation. 

In conclusion, our work substantiates the problem of the 
small mass of Mars and suggests that understanding terrestrial 
planet formation requires a paradigm shift in our view of the 
early evolution of the solar system. 
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Run 


N 


AMD 


RMC 


^planets 


^Mars 


^Mars 


Normal Migration 














Test3 1 


5 


0.0085 


32.69 


2.10 


0.56 


1.52 


Test43 


3 


0.0011 


45.51 


1.85 


0.59 


1.41 


Rani 


3 


0.0035 


58.77 


1.82 


0.43 


1.42 


Ran3 


4 


U.UUi / 


io.yj 




U.jz 


1 ^0 
l.JZ 


Ran4 


3 


0.0023 


42.09 


1.71 


0.14 


1.71 


Kan / 


J 


U.UUI 1 


jU. jU 


1. /O 


U. /o 


i.zy 


Double Asteroids 














Test61 


3 


0.0023 


42.38 


1.82 


0.59 


1.43 


Test62 


3 


0.0014 


66.43 


1.83 






Double Embryos 
Test54 


4 


0.0011 


40.53 


1.90 


0.06 


1.88 


Perfect Migration 
PM21 


4 


0.0016 


40.87 


1.93 


0.46 


1.41 


PM22 


3 


0.0113 


51.34 


1.72 


0.33 


1.37 


PM24 


6 


0.0076 


25.93 


2.06 


0.36 


1.43 


MVEM 


4 


0.0014 


89.9 


1.88 


0.11 


1.52 



Table 1. Simulation results for each simulation included are the 
number of planets Af, the angular momentum deficit (AMD) and 
the radial mass concentration (RMC) for each system of planets, 
the total mass of the planets M p i ane ts, and Mars analog MMars 
in Earth Masses, and the semimajor axis of the Mars analogs 
flMars- The entry MVEM is data for the current terrestrial planets. 
Note that Test3 1 and PM24 each had one embryo stranded in the 
asteroid belt region, with a > 2.0 AU, which are counted in their 
N, but not used to calculate RMC or AMD. 
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